Stability of spintronic devices based on quantum ring networks 
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Transport properties in mesoscopic networks are investigated, where the strength of the (Rashba- 
type) spin-orbit coupling is assumed to be tuned with external gate voltages. We analyze in detail 
to what extent the ideal behavior and functionality of some promising network-based devices are 
modified by random (spin-dependent) scattering events and by thermal fluctuations. It is found 
that although the functionality of these devices is obviously based on the quantum coherence of 
the transmitted electrons, there is a certain stability: moderate level of errors can be tolerated. 
For mesoscopic networks made of typical semiconductor materials, even cryogenic temperatures can 
smear out the desired transport properties. When the energy distribution of the input carriers is 
narrow enough, it turns out that the devices can operate close to their ideal limits even at relative 
high temperature. As an example, we present results for two different networks: one that realizes 
a Stern-Gerlach device and another that simulates a spin quantum walker. Finally we propose a 
simple network that can act as a narrow band energy filter even in the presence of random scatterers. 

PACS numbers: 85.35.Ds, 85.75.-d, 72.25.Dc 



I. INTRODUCTION 

Mesoscopic networks with Rashba-type^ spin-orbit 
coupling have recently gained considerable attention due 
to their remarkable spin transformation properties. The- 
oretical proposals for using these networks^i^iii 



specific spintronic devices often focus on the ballistic 
regime at low (zero) temperatures. Obviously, impor- 
tant features as well as fundamental aspects of the spin- 
dependent transport problem can be found using this ap- 
proach, but it is also desirable to investigate the stability 
of the results in a more realistic context. In this paper 
we focus on some promising network-based spintronic de- 
vices and study how random scatterers modify the bal- 
listic transport and investigate the extent to which the 
functionality of these devices remains close to the ideal 
one. In this aspect, the case when the input electrons are 
in thermal equilibrium at a finite temperature, i.e., they 
are not monoenergetic, is also investigated. 

The concrete model systems we have chosen to 
analyze consist of quantum rings that are made of 
semiconducting materials exhibiting Rashba-type^ 
spin-orbit interactio n^°i^^d^d3 (SOI). The strength 
of the SOI that determines the spin sensitive behav- 
ior in these rings can be tuned with external gate 
voltages^^, offering possible spintronic applications 
already with two- or three-terminal single quantum 
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Systems of quantum rings can operate as multi-purpose, 
flexible spintronic devices^, and the concept of quantum 
walk (QW) that inherently requires large systems can 
also be realized in appropriate ring networks-. In this 
paper we focus on these latter proposals, investigate 
the stability of rectangular networks of quantum rings 
as well as that of the networks that can realize QW 



with the spin degree of freedom playing the role of the 
quantum coin. 

The present paper is organized as follows. In Sec. II 
we summarize the theoretical network model to be used, 
and describe the method for taking scattering events and 
finite temperatures into account. The geometry and pos- 
sible functionality of the networks to be investigated is 
described briefly in Sec. III. Questions related to the sta- 
bility of these arrays are investigated in detail in Sees. IV 
and V. We summarize our results in Sec. VI. 



II. MODEL 

First we consider a single ring^"* of radius a in the x — y 
plane and assume a tunable electric field in the z direction 
controlling the strength of the spin-orbit interaction char- 
acterized by the parameter a^^. The Hamiltonian^^i^ in 
the presence of spin-orbit interaction for a charged par- 
ticle of effective mass m* is given by 



H = hn 



d U! 

i 1 ( 

dip 2n 



(1) 



where ip is the azimuthal angle of a point on the 
ring and the radial spin operator is given by ar/2 — 
{(TxCosip + (Ty smip)/2. The frequency associated with 
the spin-orbit interaction is denoted by w = a/ ha, while 
hO, — Ilm*cp' gives the kinetic energy with to* be- 
ing the effective mass of the electron. The position 
dependent eigenspinors and the corresponding energies 
of this Hamiltonian can be determined analytically, see 
e.g. Refs. f36||37{ for details. 

When there are leads attached to the ring that connect 
it to electric contacts, the usual question is the conduc- 
tance as a function of the SOI strength and/or the energy 
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of the incoming electrons. At zero temperature this lat- 
ter has zero variance, it equals the Fermi energy Ep of 
the system. Therefore energy conservation requires find- 
ing eigenevalues of the Hamiltonian ^ being equal to 
Ep. As a consequence of the two possible eigenspinor 
orientations as well as the possibility of the current flow- 
ing clockwise and counterclockwise along the ring, Ep 
is usually fourfold degenerate. In order to complete the 
solution of the scattering problem at a given energy, the 
spinor valued wave functions have to be joined together. 
The spinor components have to be continuous, and the 
spin current^ that enters a junction, should also leave 
it, i.e., the net current density at a certain junction has 
to vanish. Let us note that besides Griffith's boundary 
conditions'^'* that have been described above (and will be 
used throughout this paper) there are other physically 
realistic and often used possibilities as well. The choice 
of the boundary conditions is in fact shown to be related 
to the reduction of a two dimensional problem to one 
dimension^. 

The procedure above is completely analytic, the trans- 
mission matrix T that connects the input and output 
spinors as well as the conductance using the Landauer- 
Biittiker formula— can be given in closed forms. This 
method can be transferred to systems consisting of sev- 
eral rings in a straightforward way. A general scheme 
for networks to be investigated can be mathematically 
represented as 



[ l^out) 



TN\lpi, 



(2) 



thus the input state is transformed into N output spinors 
(and a refiected one, not shown above). In general the 
transmission matrices Ti , . . . , TJv are different for differ- 
ent outputs^ (see Figs. [T] and [5] for concrete geometries). 
Considering the size of these networks, it turns out to 
be practical to solve the related systems of linear equa- 
tions numerically. In this way one can determine the 
spin transformation properties of networks of consider- 
able size, and may find geometries with possible spin- 
tronic applications. 

However, in a large, or even mesoscopic system the 
transport process will not be necessarily ballistic, the 
quantum mechanical coherence of the carrier wave func- 
tions may not be maintained all along the device. In 
order to give account for this issue, we introduce inde- 
pendent point-like scattering centers in the rings. That 
is, we assume the presence of an additional potential 



(3) 



with uniformly distributed independent random posi- 
tions (f n of the Dirac delta peaks. For the sake of simplic- 
ity, the density of the scatterers is chosen such that there 



is always a single one between two neighboring junctions, 
that is, e.g. for a four terminal ring we have four scatter- 
ers. The strengths of the potentials, /7„(Z)), are random, 
they are drawn from independent normal distributions, 
with zero mean and root-mean-square deviation D. That 
is, the probability for Un{D) to have a value in a small 
interval around u is given by p{u)du, where 



p{u) = e-Kp{-u'^/2D^)/DV2n. 



(4) 



In this way, by tuning D we can model weak disturbances 
(small D) as well as the case when frequent scattering 
events completely change the character of the transport 
process (corresponding to large values of D). 

As the interference phenomena that are responsible for 
the desired behaviors of the devices are spin sensitive, we 
also introduce spin-dependent (SD) scatterers 



(5) 



where Un{D) now represents a 2 x 2 diagonal matrix, 
with independent random diagonal elements Uni{D) and 

Un2{D). 

In both cases, we use the same boundary conditions 
as at the junctions, thus we require the continuity of the 
spinor valued wave functions ^'^(i^), r{ip) at the left 
and right side of the scatterer situated at Lpn. Then the 
requirement of vanishing spin current density reads: 
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(6) 

These additional equations are still linear ones, thus 
we can find analytical solutions for any random sets 
{UniD)}, {pn}- That is, for any input state we can cal- 
culate the output spinor, that we write symbolically as 
I'^outiUniD), Pn}) ■ When after computational runs, 
the estimated density operator 

Pout[D) = 

1 (7) 

^ \^out{Un{D), Pn}) {^out{Un{D), Pn}\ , 

converges, we have all the possible information needed to 
describe what effects result from the disturbances char- 
acterized by the variable D. (Clearly, the sum above 
runs over different positions {</?„} and scatterer strengths 
{Un\, so that the independent random variables in the 
latter case correspond to normal distributions with the 
same root-mean-square deviation D.) Note that Pout{D) 
is not normalized, we can consider it as a conditional den- 
sity operator that describes the state of the electron if it 
is transmitted at all. Using the transmission probability 
T, wehaveTr[po„t(i^)]/r = l. 

In the description of the process above it was implicitly 
assumed that the incoming electrons are monoenergetic, 
but usually, at finite temperatures, this is not the case. 
Therefore, if we would like to take the energy distribu- 
tion of the input electrons also into account, we have to 



average over all possible input energies. In thermal equi- 
librium at temperature T, the output density operator 
can be written as 

Pout{T) = J p{E, T) \^out{E)) {^out{E)\ dE, (8) 

where p(£;,r) = -^[l + eii^ {E ~ Ep) /ksT]-^ . (Note 
that this expression corresponds to the Landauer- 
Biittikcr formula for the conductance at finite tempera- 
ture and low bias^.) In practice, we convert the integral 
([5]) to a sum over the possible energies, thus the expres- 
sion for Pout{T) is similar to Eq. but the weights of 
the projectors are not uniform, they are determined by 
the Fermi distribution. 
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FIG. 1: The geometry of a rectangular ring network. The 
arrows indicate the possible directions of the currents. 



III. NETWORKS AND THEIR POSSIBLE 
APPLICATIONS 

In this section we summarize briefly the properties of 
the networks the stability of which will be analyzed later. 
The first network can be seen in Fig. [U it is a rectangu- 
lar array of iV x iV quantum rings, with one input lead 
and N outputs. This geometry has already been real- 
ized experimentally^, and the conductance properties 
have also been investigated in theoretical works^. It 
has been shown in Ref. that if the strength of the SOI 
can be modulated locally (ring by ring), then relatively 
small [N = 3, = 5) networks are remarkably versa- 
tile from the viewpoint of spin transformations. Working 
in a given network geometry, the input current can be 
directed to any of the output ports, simply by chang- 
ing the SOI strengths. This kind of operation can al- 
ready be achieved by a 3 x 3 geometry in a practically 
refiectionless way, and the probability for an electron to 
leave the device through a lead other than the distin- 
guished one is less than 1%. A slightly larger network 
[N = 5) can also be made completely analogous to the 
Stern-Gerlach device: If the input is one of the eigen- 
states of (e.g. spin-up in this direction) the output 
will have the same spin direction at a certain output port 
(with probability higher than 99%). When spin direction 
of the input is the opposite (e.g. spin-down), it is directed 
towards a different output port, with its final spin direc- 
tion being the same as the initial one^. 

Having chosen a certain N x N geometry, one may 
want to realize a given spin transformation (function of 
the device). In order to achieve the desired way of op- 
eration, there are N'^ parameters to adjust, namely the 
SOI strengths in the rings. At this point, it is reason- 
able to define a real valued (preferably experimentally 
measurable) function over this iV^ dimensional param- 
eter space that has (a global) minimum when the de- 
sired goal is achieved. Using this function, an appropri- 
ate multidimensional minimum searching algorithm can 
find the relevant parameter values. (Note that in prin- 
ciple this procedure can also be applied in experimental 



situations.) Clearly, the key point here is the choice of 
the function to be minimized. It was found that rather 
straightforward definitions work properly with a "down- 
hill simplex" multidimensional minimizing routine^. If 
the purpose is to direct the input current into output 
port k, irrespectively from the spin state, then 

/i =2-Pfe(in: t)-Pfe(in: i) (9) 

is an appropriate choice for the function to be minimized. 
Here Pfc(in : f) [PA:(in : |)] is the probability of transmis- 
sion into the distinguished output port, provided the in- 
put was spin up (down) in the z direction. That is, fi is 
zero if a spin up input wave leaves the network through 
the distinguished output with 100% probability and the 
same holds for the orthogonal (spin-down) input. (Note 
that considering the two different inputs, two computa- 
tional/experimental runs are needed to obtain /i.) It is 
also worth noting that /i depends on the SOI strengths in 
a complex way, via the transmission probabilities, which 
can be measured. 

When our aim is a device that can be the spintronic 
analogue of the Stern-Gerlach apparatus, we can use the 
function 

/2 = 2 - Pi(in : T,out : T) - P5(in : i,out : i), (10) 

where Pi(in : tjOut: t) [-^5(111 : i,out: |)] denotes 
the transmission probability of a spin up (down) state 
through output port 1 (5), provided the input was spin 
up (down) . Note that the theoretical minimum of both /i 
and /2 is zero and values below 0.05 can be reached with 
SOI strengths in the experimentally achievable range. 

The second network geometry we will analyze is shown 
in Fig.O The aim of this device is the realization of quan- 
tum walk (QW), so that spin degree of freedom plays the 
role of the quantum coin, i.e., it determines the direction 
of the steps taken by the walker. The scheme itself is 
similar to the classical random walk on the line, when 
the walker tosses a coin and takes a step to the left if the 
coin is heads or to the right if the coin is tails (or vice 
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states \L) and \B) for B-type rings so that \L) {\R)) has 
zero probabiHty to be transmitted into the right (left) 
output arm. C-type rings have double role in the scheme, 
first they undo the rotation performed by the B-type 
rings (B-type rings direct the input states \L) and \R) 
into their appropriate output, but meanwhile they also 
rotate them) as well as allowing for the interference of 
subsequent steps by accepting two inputs. Using these 
three types of rings, coined QW on a line can be imple- 
mented. The actual direction of the walk is horizontal in 
Fig. [21 while the number of steps increases vertically. A 
network for QW of M steps has M -\- \ outputs and it 
consists of M(3M -I- 5)/2 elementary rings. 

Let us note that the possible QW realization summa- 
rized above uses reflectionless rings as building blocks. 
The scheme does not involve optimization, it is based on 
analytic results obtained for single rings^i^i^. All types 
of rings can be realized with several parameter pairs (size 
and SOI strength), although these optimal parameters 
form discrete sets. Therefore the system can tolerate ge- 



ometrical errors to some extent, but the optimization of 
FIG. 2: Schematic representation of a device that can realize a certain, not necessarily perfect geometry can be more 
coined quantum walk on the line. The example shown in difficult in this case than for rectangular networks, 
the figure corresponds to a QW with M = 3 steps. Note 
that in the ideal case the current flows in the network in 
the downward direction, and the input current is distributed 

among the numbered output ports. IV. STABILITY: RANDOM SCATTERERS 



versa). In the quantum analogue of such a walk a quan- 
tum coin (electron spin, in our case) is used, the state of 
which can be a linear combination of the classical heads 
and tails, or mathematically, any state of a 'coin' Hilbert 
space 'Hc^ spanned by the two basis states \L) and 
where L (R) stand for 'left' ('right'). The positions of 
the walker also span a Hilbert space Tip = {\i) : i € Z} 
with \i) corresponding to the walker localized in posi- 
tion i. The states of the total system are in the space 
Ti. = TLc ®'Hp^ which can be identified with spinor val- 
ued wave functions in our case. The conditional step of 
the walker depends on the state of the coin, it can be 
described by the unitary operation 

S = \L) {L\®Y. \' + 1) + <^I®E I* - 1) • (11) 

i i 

The coin-toss is realized by a unitary operation C acting 
on the spin states. The QW of N steps is defined as the 
transformation [/^, where U , acting on 7i = Tic ® Hp 
is given by 

U^S-{C(g)I). (12) 

The network shown in Fig. [2] realizes this scheme in the 
following manner: first A- type rings perform Hadamard 
operation on the spin states. (In other words, the op- 
erator C in this scheme is a frequently used, so called 
balanced unitary coin, which is represented by a matrix 
in which each element is of equal magnitude.) With ap- 
propriately chosen parameters, there are two input spin 



When investigating the possible effects of scatterers, 
the most obvious one is the decrease of the conductance 
of the device. This consequence of the non-ideal trans- 
port is inevitable, uncontrollable scattering events in- 
volve backscattering as well, thus the devices themselves 
can not be made reflectionless. Therefore the decrease 
of the transmission probability (or the efficiency of the 
device) as a function of the strength of the scattering 
induced disturbances is the first question to investigate. 
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FIG. 3: (Color online) The transmission probability of a sin- 
gle quantum ring with diametrically coupled input and out- 
put leads for different strengths of the scattering events. The 
strength of the SOI is chosen such that uj/n = 1.0 (see 
Eq. ©). 
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However, for some applications the most important 
point can be whether the electron's spatial and/or spin 
state is close to the desired one, provided they are trans- 
mitted at all. Reflective losses are less crucial in this 
case, we are mainly concerned in the transmitted wave 
function, and analyze to what extent it can be considered 
to be acceptable from the viewpoint of the functionality 
of the device. This brings us to the question when does 
phase randomization (i.e., the second effect of scattering) 
become important? 

Clearly, the two issues above are not independent, as 
one can check e.g. for one dimensional propagation along 
a line, stronger dephasing means stronger reflection as 
well. In a network, however, there is an additional, less 
direct connection between (partial) reflections and effec- 
tive dephasing. Namely, when the functionality of a cer- 
tain device strongly relies on multiple coherent quantum 
mechanical interferences, reflections at internal points 
modify the wave functions that interfere at the junctions, 
leading to an additional, dephasing-hke effect. 

Fig. [3] can serve as a reference, it shows for various val- 
ues of D the transmission probability of a single quan- 
tum ring as a function of the incoming wave number k 
multiplied by the ring radius a. (The input and output 
leads are connected to this ring diametrically, and the 
strength of the SOI is moderate: ui/fl = 1.0) When the 
ratio of D (measured in energy units) and the Fermi en- 
ergy is about 0.1, the transmission probability decreases 
from unity below 60%. The contrast of the minimal and 
maximal transmission probabilities also tends to disap- 
pear when D increases. For extremely strong scattering 
effects the transmission probability becomes negligible. 
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FIG. 4: (Color online) The transmission probability of rectan- 
gular NxN quantum ring arrays as a function of the strengths 
of the scattering events. The ring sizes and the SOI strengths 
correspond to the ideal operation described in the text. 

Now let us investigate rectangular arrays. The de- 
crease of the transmission probability is shown by Fig. [3] 



for both 3x3 and 5x5 arrays. The parameters used for 
these plots are the optimal ones for the ideal case (with- 
out scatterers), that is, the intended operation of the 
smaller network is the direction of the input current to 
one of its output ports, while the 5x5 array should work 
analogously to the Stern-Gerlach device. As we can see, 
in the range of moderate values of D (which is of our pri- 
mary interest), spin-dependent (SD) scatterers cause the 
transmission probabilities to decrease faster than in the 
non-SD case. A weak size effect is also visible: the same 
density of scatterers cause slightly stronger dephasing in 
a larger network. It is important to mention here, that 
although there are considerable number of rings in these 
networks, the rate at which the transmission probabilities 
decrease as a function of -D, is not orders of magnitude 
higher than it is for a single ring. As we mentioned in 
the introductory part of this section, increasing reflection 
probability is an inevitable consequence of the scattering 
events, but their presence does not necessarily lead to a 
complete destruction of the functionality of the device. 
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FIG. 5: (Color online) The measures of functionality given 
by Eq. p3p for a rectangular 3x3 quantum ring array as a 
function of the strengths of the scattering events. The ring 
sizes and the SOI strengths correspond to the ideal operation 
described in the text. 

In the case when a 3x 3 rectangular network (see Fig.[T]) 
is intended to direct its output current, various measures 
of the "distance" between its ideal and actual operation 
can be formulated in terms of the relevant probabilities: 



Ml = 

M2 = 



Pi 



Pi 



R + P2 + P3 I -Pi 

Pi ^ Pi 
P2 + P3 T~Pi 

^lo 



(13) 



k 



where R and T denote the reflection and transmission 
probabilities, and Pk stands for the probability that the 
electron leaves the device through output arm fc (fc = 1 
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corresponds to the distinguished output.) The measures 
Ml and M2 are signal to noise ratios, with Mi being the 
more strict one in the sense that M2 does not take reflec- 
tive losses into account; M2 compares the probability of 
the desired output to that of the unwanted ones. Sout is 
the Shannon entropy of the output, it should be zero in 
the ideal case, and has a maximum when all Pk have the 
same values. The dependence of these measures on the 
strength of the disturbances caused by the scatterers is 
shown in Fig. [S] As we can see, these measures decrease 
faster than the transmission probability, but acceptable 
signal to noise ratios can be seen even when the transmis- 
sion probability is below 80%. Note that the difference 
between effects caused by SD and non-SD scatterers is 
surprisingly small, they destroy quantum mechanical co- 
herence at a similar rate. This is a general result for all 
the investigated devices, in the following we will focus on 
the spin-dependent case. 



10 



1 r 



0.01 



-M'j, SD scatterers 

■ M'^, SD scatterers 

■ S' , SD scatterers 



0.00 



D/E^ 



FIG. 6: (Color online) The measures of functionality given 
by Eq. (|14|l for a rectangular 5x5 quantum ring array as a 
function of the strengths of the scattering events. The ring 
sizes and the SOI strengths correspond to the ideal operation 
described in the text. 

Next we focus on the spintronic analogue of the Stern- 
Gerlach device proposed in Ref. p|]. According to the 
functionality of this device, the relevant measures are 
spin-dependent : 



M'l = 



M'2 



Pi 



TTi 



TTi 



TTi 



E 



Em. 

T 



(14) 



l0g2( 



TTfe- 



T 



T 



T 



where we considered the case of spin- up input, and 
P|-|-i = Pi (in : t, out : t) is the probability of the desired 



TT 

output. Pi 



Shannon entropy is related to spin-down output at the 
channel k for spin- up input. (Note that similar measures 
can be defined also for the orthogonal input, but due to 
symmetry, those measures give the same numerical values 
as the ones above.) 

Fig. [6] shows the dependence of the measures given by 
Eq. ([13]) on D. Note that the limiting value of S'^^^ for 
strong scattering events (and also that of Sout shown in 
Fig. [5]) is essentially the maximum of the Shannon en- 
tropy, i.e., in this limit the (weak) output current is dis- 
tributed evenly among the output channels. Comparing 
Figs. inland [5] it can be seen that the same value of D in- 
duces considerably stronger effects in the larger network, 
although even the latter can be functional for a moderate 
level of disturbances. Size effect, however, can be ana- 
lyzed best in QW networks, where the desired function- 
ality is the same, it is only the number of the constituent 
rings that is different for networks that are designed to 
realize different number of QW steps. 

In the case of the quantum walk, however, determin- 
ing the spatial current distribution at the output (i.e, the 
probabilities corresponding to the transmission through 
the respective outputs) is not sufScient to quantify the 
effectiveness of the device: phase relations are also of fun- 
damental importance. Therefore we calculate the over- 
laps 



Pi = Tr{poutPi), F2 = 



Tr(po„fp.,) 
T 



F. = 



Tr(po„t/Jt 
T' 



Pfc(in : t,out : i) in the spin resolved 



(15) 

where pout is given by Eq. ([7]), and pi is the density op- 
erator corresponding to the pure quantum mechanical 
state that would result from the ideal QW described by 
Eq. In contrast to the transmission probability T, 

that does not distinguish the outputs, T' denotes the 
probability of transmission into "valid" output channels: 
In the ideal case there is no output into the leads at 
the sides of the device, the input current is distributed 
among the output leads located at the bottom side of 
Fig. [51 However, scattering effects modify this picture, 
and T — T' , the transmission probability into the side 
leads (that is zero in the ideal case), can also be con- 
sidered as a measure of the functionality of the device. 

Considering transmission probabilities T of networks 
that are designed to realize QW, the most remarkable 
point is that the functions T{D) are quite similar for 
different network sizes (relative errors are below 5% for 
M = 1,...,9 QW steps). If we recall that the input 
part of these devices of any size is exactly the same, the 
conclusion is that a small subnetwork close to the input 
junction is responsible for the reflection, i.e., scattering 
events in this small area determine the reflection proba- 
bility P = 1 — T. The measures Pi, P2 and P3, however, 
show stronger size dependence, thus the current is dis- 
tributed among the output channels in different ways for 
networks of different sizes. The overlaps given by Eq. (jlSp 
are shown in Fig. [71 as a function of the strength of the 
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FIG. 7: (Color online) The measures of functionality given 
by Eq. psp for quantum ring arrays to realize A/ step QW. 
The ring sizes and the SOI strengths correspond to the ideal 
operation described in the text. 



scattering events. Recalling that the number of the el- 
ementary rings in a network to realize M QW steps is 
proportional to M^, one might expect the overlaps Fk 
measuring the fault tolerance of the devices to scale also 
as M^. However, choosing a certain level of any of these 
measures (e.g. F2 — 0.9), we find that the D values cor- 
responding to this level do not decrease quadratically as 
the size of the networks increases. Instead, this decrease 
is even weaker than a linear dependence, thus the sta- 
bility of the devices scale with size in a promising way. 
Note that similar error analysis for the possible quan- 
tum optical realizations of QW has been investigated in 
Refs. giE^. 



V. STABILITY: TEMPERATURE RELATED 
EFFECTS 

The role of finite temperatures are taken into account 
here as the appropriate broadening of the energy distri- 
bution of the incoming electrons. In view of this, there 
is a natural limit, describing the case when interference 
related effects are expected to appear at all. Namely, if 
there is a typical wave number k, and a characteristic 
length a in a quantum interference device, then interfer- 
ence properties are (quasi)periodic as a function of ka. 
(In the case of a single ring with moderate SOI strength, 
k can be the wave number of the input wave, which is 
assumed to be monoenergetic at this point, while a is 
the radius of the ring.) Obviously, if the incoming wave 
numbers have a range of Afc, aAk cannot be too large 
(i.e., comparable to unity), otherwise the interferences 
are smeared out. As we can see, there is a size effect 
here: at a given temperature T, the width of the equi- 
librium energy distribution is fixed, it has the order of 
ksT, and Afc is also given by the dispersion relation. 
Thus if we can decrease aAk e.g. by building smaller de- 



vices, then temperature related effects will become less 
pronounced. In other words, miniaturization, besides its 
obvious advantages, can help to avoid interferences to 
disappear (and also decrease dephasing caused by ran- 
dom scattering, simply by increasing the mean free path 
- system size ratio). 




20.0 20.5 21.0 21.5 22.0 

ka 



FIG. 8: (Color online) The transmission probability of a single 
quantum ring with diametrically coupled input and output 
leads for different temperatures. The strength of the SOI is 
characterized by u/Q = 1.0 (see Eq. (H}). 

Before we analyze more complex systems, let us con- 
sider again a single quantum ring with one input and one 
output arm, with ka being in the experimentally achiev- 
able range around 20. Fig. [8] shows ~ in accordance with 
the introductory part of this section - that the transmis- 
sion probabilities are noticeably modified already when 
ksT is a small portion of the Fermi energy. (Numer- 
ically, for Ep being in the range of 10 meV, the flat 
line that belongs to the highest temperature in Fig. [51 
is still in the cryogenic range.) Comparing the two lim- 
its when non-monoenergetic input and strong scattering 
events completely change the transmission probabilities 
(see Figs. [5] and [3]), an important difference becomes 
apparent: high temperatures lead to finite transmission 
probabilities, while strong scattering events cause the de- 
vice to be opaque. The latter result can be understood 
readily, but the first one may need some explanation. 
High temperature conductance is an integral property in 
the sense of Eq. ([8]), if the input is normalized, we may 
write^ 

G(T) = ^ j p{E,T){-^out{E)\^out{E))dE. (16) 

Thus all monoenergetic transmission probabilities con- 
tribute to the conductance with a weight determined by 
the energy distribution. Recalling that the transmission 
probability is a quasiperiodic function of ka, for wide 
enough input energy distribution (being almost constant 
during a period) the averaged conductance is propor- 
tional to the integral of the transmission probability over 
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a period divided by the length of the period. For a given 
geometry, this quantity is a constant, but it can change 
when the strength of the SOI is modified. 

Fig. [5] shows that a 3 x 3 rectangular array can direct 
its input current to one of its output ports with M2 — 2 
signal to noise ratio at a temperature, where the trans- 
mission probability of a single ring reduces already below 
80% of its initial value (see Fig. [S] as well.) The spintronic 
analogue of the Stern- Gerlach device based on a 5 x 5 
network is more sensitive to thermal fluctuations, which 
is due to its larger size and different functionality. (It 
is worth noting that the Shannon entropies reach their 
maxima in the high temperature limit as well.) Accord- 
ing to Fig. [ini the stability of the QW networks scale 
with their size M similarly to the case of random scat- 
terers: given values of F2 is reached at temperatures that 
decrease as a function of M slower than linearly. 



10 - 



M\ (N=5) 



■M^ iN=3) 
-S„,., iN=3) 



SD scatterers 



0.00 




0.04 



FIG. 9: (Color online) The measures of functionality given by 
Eqs. (|13|) and (|14|l for rectangular 3x3 and 5x5 quantum 
ring arrays as a function of temperature. The ring sizes and 
the SOI strengths correspond to the ideal operation described 
in the text. 
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FIG. 10: (Color online) The measures of functionality given 
by Eq. I|15p for quantum ring arrays to realize M step QW. 
The ring sizes and the SOI strengths correspond to the ideal 
operation described in the text. 



be the most effective method to overcome the problem 
caused by thermal broadening of the input energies, but 
it can serve demonstrative purposes.) We found that 
an ideal linear chain consisting of five rings of linearly 
increasing sizes (-1-0.5% for the consecutive rings) has 
a transmission peak the width of which corresponds to 
ksT = 10^^ X Ep. The transmission probability is es- 
sentially zero in an energy interval of ks/Ep — 0.06 (At 
the peak it is around 60%.) These ideal properties are of 
course modified with introducing scatterers like we did 
in the current paper. However, it was found that for 
moderate strength of the scattering effects (e.g. when the 
overlap F3 is around 0.5 for a nine step QW), the height 
of the transmission peak is still above 10%, and its char- 
acteristic width corresponds to fc^T = 10~* x Ep, which 
is remarkable as we did not optimize the geometry for 
these values. 



These results, however, still mean that quantum inter- 
ference devices require quasi-monoenergetic input elec- 
trons, otherwise interference effects are smeared out. In 
usual semiconductors, the limit temperature is quite low, 
meaning an obstacle in many practical applications. This 
effect can be decreased by building smaller devices, or 
decreasing the width of the input energy distribution be- 
low the thermal equilibrium. Note that in this latter case 
nonelastic scattering events will eventually lead to ther- 
mal equilibrium again, but the thermalization length be- 
ing relevant here can be considerably larger than e.g. the 
dephasing length. 

As a final point, we investigated whether simple 
network-based devices can provide sufficiently narrow 
transmission peaks, i.e., whether they can deliver quasi- 
monoenergetic electrons. (We note that this might not 



VI. SUMMARY 

In this paper we investigated the stability of network- 
based spintronic devices against random scattering 
events and thermal fluctuations. We focused on the ques- 
tion to what extent the functionality of certain proposed 
networks remains close to the ideal behavior when tem- 
perature or the intensity of the scattering induced dis- 
turbance increases. We considered two different network 
types, a rectangular array that can serve as a spintronic 
analogue of the Stern-Gerlach apparatus as well as net- 
works that can realize quantum walk with the spin degree 
of freedom playing the role of the quantum coin. We 
introduced appropriate measures to quantify the "dis- 
tance" between the ideal and the realistic behavior of 
these devices, and found that the functionality can tol- 
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erate moderate level of errors. It was demonstrated that 
although networks of increasing size are less stable, but 
fortunately this decrease of stability is slow with increas- 
ing the size of the networks. 
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